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Abstract 

Highly excited many-particle states in quantum systems (nuclei, atoms, quan- 
tum dots, spin systems, quantum computers) can be "chaotic" superpositions 
of mean-field basis states (Slater determinants, products of spin or qubit 
states). This is a result of the very high energy level density of many-body 
states which can be easily mixed by a residual interaction between particles. 
We consider the time dynamics of wave functions and increase of entropy in 
such chaotic systems. 

As an example we present the time evolution in a closed quantum com- 
puter. A time scale for the entropy S{t) increase is tc ~ ro/(nlog2n), where 
To is the qubit "lifetime", n is the number of qubits, 5(0) = and S{tc) = 1. 
At t <^tc the entropy is small: S ~ nt^ log2(l/t^ J^), where J is the inter- 
qubit interaction strength. At t > tc the number of "wrong" states increases 
exponentially as 2'^^*^ . Therefore, tc may be interpreted as a maximal time for 
operation of a quantum computer, since at t > tc one has to struggle against 
the second law of thermodynamics. At t 3> tc the system entropy approaches 
that for chaotic eigenstates. 
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I. INTRODUCTION 



Highly excited many-particle states in many-body systems can be presented as "chaotic" 
superpositions of shell-model basis states - see the recent calculations for complex atoms [|l| , 
multicharged ions 0, nuclei [|| and spin systems 0,^. Indeed, the number of combinations 
to distribute n particles over m orbitals is exponentially large (m!/n!(m — n)! in a Fermi 
system). Therefore, the interval between the many-body levels D is exponentially small and 
residual interaction between the particles mixes a huge number of the mean-field basis states 
(Slater determinants) when forming eigenstates. The number of principal basis components 
in an eigenstate can be estimated as A'p ~ F/D where F is the spreading width of a typical 
component that can be estimated using the Fermi Golden Rule. In such chaotic eigenstates 
any external weak perturbation is exponentially enhanced. The enhancement factor is ~ 
\Jn~p oc 1/a/D - see e.g. and references therein. This huge enhancement have been 
observed in numerous experiments studying parity violation effects in compound nuclei - see 
e.g. review and references therein. 

In a recent work |^ the consideration of many-body chaos has been extended to quantum 
computers [p|-p^. Any model of a quantum computer is somewhat similar to that of a spin 
system. In Ref. the authors modelled a quantum computer by a random Hamiltonian, 

^ = Ee^^r + E^^X^P (1) 

where the cTj are the Pauli matrices for the qubit i and the second sum runs over nearest- 
neighbor qubit pairs. The energy spacing between the two states of a qubit was represented 
by ej which was uniformly distributed in the interval [O.SAq, I.SAq]. Here e, can be viewed as 
the splitting of nuclear spin levels in a local magnetic field, as discussed in recent experimen- 
tal proposals 0)111] • The different values of e, are needed to prepare a specific initial state 



by electromagnetic pulses in nuclear magnetic resonance. In this case the couplings Jij will 
represent the interactions between the spins, which are needed for multi-qubit operations 
in the quantum computer. The total number of states in this system is = 2^^, and the 
typical interval between the nearby energies of multiqubit states is ~ Aon2~"'. 

A rough estimate for the boundary of the chaos in the quantum computer eigenstates 
is Jc ~ Ao/qn, where qn is the number of interacting qubit pairs {qn = 2?t, in a 2D square 
array of "spins" with only short-range interactions) . This follows from a simple perturbation 
theory argument: the mixing is strong when the perturbation is larger than the minimal 
energy interval between the basis states which can be directly mixed by this perturbation 
(see detailed discussion of the boundary of chaos in many-body systems in Refs. [p3|-p5|]). 
Numerical simulations in [|| have shown that the boundary of the chaos in the quantum 
computer eigenstates is Jc — 0.4Ao/n . Above this point they observed a transition from 
Poissonian to Wigner-Dyson statistics for the intervals between the energy levels. For J < 
one eigenstate is formed by one or few basis states built from the non-interacting qubits 
(products of "up" and "down" states). For J > Jc a. huge number of basis states are 
required. 

Because of the exponential laws it is convenient to study the entropy S of the eigenstates 
(in many-body systems the entropy 5* ~ InA'p, see e.g. Ref. |2^). In Ref. they observed 



c 



a dramatic increase of the eigenstate entropy in the transition from J < Jc to J > Jc] in 
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fact , they defined Jc as a point where S = 1. In Ref. ^ this process of chaotization of 
the eigenstates with the increase of J, or number of qubits n, was termed as a melting 
of the quantum computer and was assumed to lead to destruction of its operability. The 
authors stress that this destruction of operability takes place in an isolated (closed) system 
without any external decoherence process ( one could complement this picture by the \jNp 
enhancement of any weak external perturbation acting on the quantum computer). 

This straightforward conclusion may be misleading. "Theoretically" , this picture is sim- 
ilar to that observed in nuclei and atoms. However, the "experimental" situation is very 
different. In nuclei and atoms experiments have resolved particular many-body energy lev- 
els. Therefore, the description of the systems based on a consideration of the eigenstates 
was an adequate one. In quantum computers the energy interval between the eigenstates 
is extremely small. The authors of Ref. [§] estimated that the average interval between the 
multi-qubit eigenstates for 1000 qubits, the minimum number for which Shor's algorithm 
10| becomes useful [jl3[ is D ~ 10~^^^ii' (for a realistic Aq ~ IK). Therefore, in the case of a 



quantum computer it is impossible to resolve multiqubit energy levels. Temperature, or the 
finite time of the process r, gives an uncertainty in energy 5E ^ D. In this case the picture 
with chaotic eigenstates is not an adequate one and we should consider the time evolution 
of the quantum computer wave function and entropy. Quantum chaos in the eigenstates 
allows us to apply a statistical approach to this consideration. 



II. TIME EVOLUTION OF THE CHAOTIC MANY-BODY STATE 

Exact ("compound") eigenstates |A;)of the Hamiltonian H can be expressed in terms 
of simple shell-model basis states |/) in many-body systems or products of qubits in a 
computer: 

l^) = Ecfl/); |/) = a+...a;jO). (2) 

/ 

These compound eigenstates \k) are formed by the residual interaction J ; are creation 
or spin-raising operators (if the ground state |0) corresponds to spins down). Consider now 
the time evolution of the system. Assume that initially (t = 0) the system is in a basis state 
\i) (quantum computer in a state with certain spins "up") which can be presented as a sum 
over exact eigenstates: 

K) = Ecf^|fc)- (3) 

k 

Then the time- dependent wave function is equal to 

^W = EcfC'f |/)exp(-.E(^)t). (4) 

The sum is taken over the eigenstates k and basis states / ; we put % = 1. The probability 
Wi = l^jp = I (i|\E'(t)) p to find the initial state in this wave function is determined by the 
amplitude 
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Ai = {i\ exp{-iHt)\i) = J2 |Cf ^pexp(-zE('=)t) ~ / dEPi{E) exp{-iEt). (5) 

k 

Here we replaced the summation over a very large number of the eigenstates by the inte- 
gration over their energies E = E^^"* and introduced the "strength function" Pi{E) which is 
also known in the literature as the "local spectral density of states" , 



pm^\cr'?p{E), (6) 

where p{E) is the density of the eigenstates. In chaotic systems the strength function is 
given by a Breit-Wigner- type formula PB| , P5| : 



p.(E) = — (7) 



r,{E) ~ 2n\H,f\'pf{E) ~ JV/Ao- (8) 

Here 6i is the correction to the unperturbed energy level Ei due to the residual interaction J, 
Pf{E) ~ qn/Ao is the density of the "final" basis states directly connected by the interaction 
matrix element Hif with the initial state We see from the equations above that the time 
dynamics is determined by the structure of the eigenstates. 

It is easy to find Wi{t) for a small time t. Let us separate the energy of the initial state 
Ei = Ha in the exponent and make a second order expansion in H — Ei ot E — Ei in eq.(§). 
The result is: 

Ai = exp{-iE,t) (1 - {AE)Hy2), (9) 
W,{t) = (1 - {AE)H'), (10) 

{AEr = T.Hfj = T.4 = 1^Jr- (11) 



Here (AE)"^ is the second moment of the strength function, is the r.m.s. value of the 
interaction strength, = Jfy The first moment is equal to Ei = Ha, see e.g. |^J2 



\2 

IS LliC SCUI 

t2 = "a ^^^^ ^^^^^ ^^^^^ 

where one can also find the calculations of (AE)"^ and spreading width T{E) for many-body 
systems. 

Note that in the special case of a very strong residual interaction J ^ Aq this short-time 
dependence can be extended to a longer time using the exact solution for the case of Aq = 
(in this case it is easy to calculate exp{—iHt)): 

Ai = n cos J,,t - 11(1 - iJ^Jtf/2) exp{-{AE)hy2), (12) 

i<j i<j 

W^,(t) ~exp(-(AE)V). (13) 
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The strength function and density of states in this hmit are also described by Gaussian 
functions with variance cr^ = (AE)"^: 



2" 

'"^> = 72iP^^'"-^>- '''' 

The density of states remains Gaussian for Aq 7^ 0, with cr^ = ne^ + [AE^ if there is no gap 
in the single-qubit spectra (in Ref. 0] the "up" and "down' spectra were separated by a gap 
equal to Aq). In general the unperturbed density of states (J = 0) can be presented as a 
sum of the Gaussian functions (one should separate classes of states with a certain number 
of spins "up"). The interaction J in the Hamiltonian (|l|) mixes these classes and makes the 
density closer to the single Gaussian function. 

The limit at large time in chaotic case can be obtained by calculation of the integral in eq. 
(^) in the complex E plain. We should close the contour of integration in the bottom part of 
the complex plane {Im{E) < 0) to provide a vanishing contribution at infinity. The limit at 
large time t is given by the pole of the strength function (|^ closest to the real E axis . If F 
and 6i do not depend on E the integration gives the usual exponential decay Wi = exp(— Ft) 
26| . However, the dependence of the spreading width on energy E is necessary to provide 



the finite second moment (AE)"^ of the strength function (note that in many-body systems 
the dependence T{E) can be approximated by a Gaussian function, since the density of final 
states Pf{Ef) in eq. is usually close to Gaussian [0). If F < AE the closest pole is given 
by F = —2Im{Ep), where Ep is a solution of the equation Ep = Ei + 6i{Ep) — iT{Ep)/2 with 
a minimal imaginary part. If F -C AE we have F = F. As a result we obtain an exponential 
dependence for large t: 

lyi(t) ~exp(-n). (16) 

It is useful to have a simple extrapolation formula (valid for F < AE) between the cases of 
small time eq. (0) and large time eq. (|T6|) : 



Wi{t) = exp 



F 



2 



2(AE)2 \ A{AE) 



+ TH^ I . (17) 



Now we can estimate the probabilities of the other components Wf. For small time or 
small interaction J, other components can be populated due to direct transitions from the 
initial state only: 



Wf = \{f\exp{-iHt)\i)\'^\H,f\ 



|2 



Ai{t) \ exp{iiUift)dt 
exp {iuJif — F/2)t-l|^ (18) 
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Here ooij = Ej — Ei. We stress again that this approximate equation does not contain 
transitions between the small components. For example, it does not contain the width of 
the state /; the width F stands only to indicate some increase of the denominator and to 
clarify the "small time" condition that should include small Uift, Ft or AEt. For small 
time Wf = \Hif\'^t'^. Here Hij is equal to one of the Jij that produces a change of the state 
of a pair of "spins" (qubits), transferring initial state i to another state /. The result at 
larger times is different for perturbative and chaotic regimes. In the perturbative regime, 
J -C Ao/qn eq. ([T8|) is the final one. In the chaotic regime we can find the asymptotic 
expression for large times. The projection of in eq. (^) to the component / gives 

Wf{t) = W' + WP'''\t), (19) 



Here F* ~ Fi + F/ ~ 2F. 

WP'''\t)= J2 Cf^Cf Cp^Cf exp(2(E(^)-£;(p))t). (21) 



At large time t, the different terms in Wj ^'^ (t) rapidly oscillate and we can put Wf^^'^t) = 
0. Thus, asymptotically the distribution of the components in the time-dependent wave 
function is close to that in the chaotic eigenstates (see eqs (|^,|^)) with a doubled spreading 
width. 



III. ENTROPY INCREASE 

It is convenient to define the entropy of a many-body state as a sum over the basis 



components ( a comparison with other definitions can be found, e.g. in Ref. p^): 



S = -J2Ws log2 Ws = -W, log2 W,-^Wf log2 Wf. (22) 
Initially, we have only one component, Wi = 1 and the entropy is equal to zero. It is easy 



to obtain a small-time estimate for the entropy using eqs. (p!0|, [IT], IT 



S ^ {AEYt' log,{qn/{AE)V) = qnJ^t' log,{l/ J^t'). (23) 

We see that the initial increase of the entropy is relatively small (~ t^), however, it is 
proportional to the number of qubits n. 

The criterion of a quantum computer "melting" used in Ref. is the entropy 5 = 1. We 
can extend the small- time consideration to include this point. For small time we have some 
decrease of the initial component and population of the components directly coupled to the 
initial one. The number of such small components is equal to the number of interacting 
pairs (qn) in the Hamiltonian eq.(|lD, since each pair can change its state due to interaction 
and this leads to a different many-body state. Using the normalization condition J2s = 1 
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we obtain an estimate Wf = (1 — Wi)/nf where Uf is the "principal" number of the final 
components. Initially ^/ = qn. This gives us the following approximate expression for the 
entropy: 

S = -W, log2 W, - log2((l - W,)/nf) ^ Wf 
= -W, log2 W,-{1- W,) log2((l - W,)/nf) ~ (1 - W,) log^irif). (24) 

The last approximate expression is an estimate with logarithmic accuracy, assuming log2(rij) 
is large. 

The condition S = 1 combined with eq.(0) for Wi{t) and eq.(|2^ for the entropy S{t) 
gives 



Wi{t) = exp 



y2{AEy \\ A{AE) 



+TH^\ =l-l/\og,{nf). (25) 



This means that the "melting" happens when the probability to be in the initial state Wi is 
still close to 1 (since log2(r;,/) is large). The loss of operability of the quantum computer is 
due to the admixture of a large number of the small components ("wrong" basis states). 

We should note that , strictly speaking, the argument of the log2 may differ from rij = qn, 
since the point can be outside the small time approximation. However, the estimate in eq. 
(P5|) with log2?T./ — log2n is valid with logarithmic accuracy (for example, a more accurate 
estimate in the case of F ^ AE is ~ qnT /Aq; this follows from eq.([T8D). 

Equation (|25|) allows us to obtain a simple estimate for the maximal operational time tc'. 



tr 



logo n 

(26) 



- riog2(n)^ (AEy • 

In the case of F -C AE we have 

~ ^ = I2 (27) 

riog2(n) n\og2{n) 

Here tq = h/To is the "lifetime" related to a single qubit, Fq = T/n; recall that F is 
proportional to the number of qubits n. More accurate result can be obtained numerically 
using expressions for Wi and Wf presented above. 

At this point we can say something about the effects of the environment. They also 
lead to "depolarization" of a qubit, which means nonzero probability of the opposite spin 
state. If this probability is small we can speak about the probabilities of the population of 
n many-qubit basis states. Each admixed basis state in this case has one of the qubit states 
different from the initial state. To account for this effect one may use a real (experimental) 
qubit lifetime tq in the estimate (pTf). 

For t > tc the higher orders in Hfft'^ expansion become important and the number of 
the small components increases exponentially: each state generates qn new states. This 
corresponds to an approximately linear increase in the entropy. At t ^ tc we can use the 



asymptotic form (pO]) of the component distribution. It is two times broader (F^ = 2F) 



than the basis component distribution of chaotic stationary states. This means that the 
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asymptotic number of the principal components is equal to Np(t) = 2N^\ where N^^ ~ 
r/D is the number of principal components in a chaotic eigenstate. It is easy to calculate 
the entropy in this case. From the normalization condition = 1, it follows that 

WT = l/iVp. Then 

S=-Y.Ws log2 Ws ~ log2 NpY.Ws = log2 Np. (28) 

s s 

Thus, the asymptotic value of the entropy is S{t ^ tc) = log2(2iV^'^)) = 5'*^'^^ + 1, where 
S^'^^ = log2 A^^'^-' is the entropy of a chaotic eigenstate. Note, that it is smaller than the 
maximal possible entropy Smax = log2 2" = n. This is due to localization of the wave 
function within the energy shell centered at the energy of the initial state with the width 
2T. 

IV. CONCLUSION 

The time dependence of the closed quantum computer wave function is different in the 
non-chaotic and chaotic regime. In the non-chaotic case J <^ Ao/n, the number of principal 
components A'p ~ 1 and the wave function remains localized near the initial state (as it was 
pointed out in the energy level density of the many-qubit states can be exponentially 
high even in this case). An increase in the number of qubits n leads to a transition to a 
chaotic regime where J > Aq/u. In this case one can operate the quantum computer within 
a limited time t < tc = TQ/nlog2n, where tq is the "lifetime" of one qubit. For t > tc it 
is hardly possible to operate the quantum computer, since in this case one faces a hopeless 
struggle against the second law of thermodynamics: increase of the entropy S{t) and very 
fast exponential increase of the number of "wrong" states Np{t) = exp S{t). The asymptotic 
value of the entropy is then close to that for chaotic eigenstates. 

A similar picture for the entropy increase is expected in other many-body systems. For 
example, one can consider a decay of a single-electron wave function in a many-electron 
quantum dot. In this case tc ~ r/log2n/ where ^/ is the effective number of final states 
that contribute to the decay width F = h/r. One may also speculate about the "entropy" 
increase for decay of a single-particle wave function in chaotic quantum billiard or disordered 
system using expansion of this wave function in the plane wave basis or the orbital angular 
momentum basis. 

This work was supported by the Australian Research Council. The author is grateful 
to M.Yu. Kuchiev for valuable discussion and to A.S. Dzurak for careful reading of the 
manuscript. 
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